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In this paper the improved multipolar model of bremsstrahlung accompanied the a-decay is presented. 
The angular formalism of calculations of the matrix elements, being enough complicated component of the 
model, is stated in details. A new definition of the angular (differential) probability of the photon emission 
in the a-decay is proposed where direction of motion of the Q-particle outside (with its tunneling inside 
barrier) is defined on the basis of angular distribution of its spacial wave function. In such approach, 
the model gives values of the angular probability of the photons emission in absolute scale, without its 
normalization on experimental data. Effectiveness of the proposed definition and accuracy of the spectra 
calculations of the bremsstrahlung spectra are analyzed in their comparison with experimental data for the 
^^"Po, ^^"^Po, ^■^''Ra and ^''^Cm nuclei, and for some other nuclei predictions are performed (in absolute 
scale). With a purpose to find characteristics taking infiuence on the bremsstrahlung probability strongly, 
a dependence of the bremsstrahlung probability on effective charge of the decaying system is analyzed. As 
a natural result, as supposed type of decaying system the emission of proton from nucleus is studied, for 
which the effective charge is essentially larger in a comparison with a-decay. For some proton emitters 
estimations of the bremsstrahlung probability are obtained (at first time, in fully quantum approach). Also 
the bremsstrahlung in fission of the ^^^Ci nucleus is analyzed (at first time, in fully quantum approach). 
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Introduction 



j3^or last two decades many experimental and theoretical efforts have been made to investigate on the nature of 
the bremsstrahlung emission accompanying a-decay of heavy nuclei. A key idea of such researches consists in 
finding a method of extraction of a new information about dynamics of a-decay (and dynamics of tunneling) 
from measured bremsstrahlung spectra. 

The first paper [I] at this topic was published by /. S. Batkin, I. V. Kopytin and T. A. Churakova as 
far back as in 1986 where a general quantum-mechanical formalism of the calculation of the bremsstrahlung 
spectra in the a-decay was constructed and the bremsstrahlung spectrum for ^^'^Po was estimated in the relative 
scale for photons energies up to 200 keV (until the fulfillment of the first experiments). In 1991 C. J. Luke, 
C. A. Gossett and R. Vandenbosch [2\ measured 7-rays for spontaneous fission of the ^'^^Cm nucleus and 
established upper limit for photon emission. In 1994 A. D'Arrigo, N. V. Eremin, G. Fazio, G. Giardina et al. 
[3] at first time measured bremsstrahlung emission in a-decay obtaining spectra for the ^^^Po and ^^^Ra nuclei. 
They analyzed the spectra in frameworks of instant accelerated models constructed on the basis of classical 
electrodynamics with a purpose to study dynamics of a-decay. In 1996 M. I. Dyakonov and I. V. Gornyi 
in [4] studied bremsstrahlung processes during tunneling of electromagnetic charge inside the potential barrier. 
They constructed a main basis of semiclassical approach for calculation of the spectra of bremsstrahlung during 
a-decay. In 1997 J. Kasagi, H. Yamazaki, N. Kasajima, T. Ohtsuki and H. Yuki in [51 [S] published the first 
measurements of the bremsstrahlung accompanying a-decay of the ^^'^Po and ^^^Cm nuclei using Ge detectors. 
They analyzed the spectra on the basis of calculations of the emission probability in the frameworks of model 
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of M. I. Dyakonov and I. V. Gornyi. They found that the bremsstrahlung probabihty is mach sniaUer then 
estimated by Coufomb accelerated modeL A natural and important conclusion from such estimation was the 
amplitude of emission of photons from the barrier cannot be neglected, that opens importance of direct quantum 
study of photon emission during tunneling. 

In 1998 T. Papenbrock and G. F. Bertsch in [7J constructed directly quantum model of calculation of 
the bremsstrahlung spectra in the a-decay developed on the basis of quantum electrodynamics with use of 
perturbation theory. Here, wave function (current) of photons is described in the dipole approximation, matrix 
element is calculated with application of the Fermi golden rule. They estimated the bremsstrahlung spectrum 
for ^'^'^Po. This paper essentially improves the previous formalism in X and its important achievement is in a 
possibility to calculate bremsstrahlung spectra with enough good accuracy. And further it was often used as a 
main basis by other independent research groups. 

In 1999 A'^. Takigawa, Y. Nozawa, K. Hagino, A. Ono and D. M. Brink [8^ performed quantum analysis of 
the bremsstrahlung in the a-decay of ^^"Po on the basis of semiclassical approach. Introducing idea of mixed 
region, they studied the contribution from it, from tunneling, outside barrier regions and from the wall of the 
inner potential well to the total spectrum. They analyzed a possibility to eliminate the ambiguity in the nuclear 
potential between the a-particle and daughter nucleus using the bremsstrahlung spectrum. Practically, at the 
same time but independently E. V. Tkalya [9l [TO] constructed the multipolar approach in frameworks of the 
fully quantum model proposing three different ways of calculation of the bremsstrahlung spectra. Such approach 
opens a possibility to study in this problem the next corrections of bremsstrahlung after dipole one being more 
accurate in a comparison with the dipole quantum approach of T. Papenbrock and G. F. Bertsch. Calculating 
spectra for ^^°Po, ^^''Po and ^^^Ra, he at first time estimated the multipole correction E2. C. A. Bertulani, 
D. T. de Paula and V. G. Zelevinsky in |11] studied tunneling of a-particle by directly solving the non- 
stationary Schrodinger equation and on such basis obtained the bremsstrahlung spectrum. They found the 
large deviation from the classical bremsstrahlung spectrum. M. I. Dyakonov in T2] found a classical formula of 
the bremsstrahlung spectrum in a-decay and showed that taking tunneling into account is crucial in accurate 
description of the bremsstrahlung. In 2000 W. So and Y. Kim, calculating the probability of the photons 
emitted in a-decay of the ^^"Po, ^^''Po nuclei by the approach of T. Papenbrock and G. F. Bertsch and at first 
time predicting the spectra for the ^^^Po, ^^"Rn, ^"'"Cm, 248(^f nuclei (in absolute scale), the first studied a 
dependence of the bremsstrahlung probability on the a-particle energy and the daughter nuclear charge |13j . 
Such analysis allows to consider deeper the a-decay and can be useful in search of nuclei with the maximal 
bremsstrahlung that can be useful in future experimental study of bremsstrahlung in a-decay. 

In [14|, I15j discussions between both experimental groups concerning measurements of the bremsstrahlung 
spectrum for the ^^'^Po nucleus were published. Here, a main problem was in the following: 1) spectrum obtained 
for ^^°Po in [5] is smaller essentially in a comparison with spectra obtained for ^^^Ra and ^^''Po in 3^; 2) unclear 
difference between the spectra trends is observed: while spectra for ^^^Ra and ^^"'Po in [3 are monotonous, 
spectrum for ^'^^Po in [5] has non-monotonous trend with a "hole" at 275 keV. The authors found explanation 
lied in the different values of angle between photon emission and a-particle motion used in experiments. Such 
discussions had caused an increased interest. If the difference between the bremsstrahlung angular spectra will 
be confirmed then the bremsstrahlung intensity should depend on the directions of the photons emission and 
the a-particle motion. So, a three-dimensional picture of the a-decay with the bremsstrahlung in the spatial 
region of nuclear boundaries has been devised. However, if the bremsstrahlung varies visibly with changing of 
the angle value, then one can suppose that the photon emission is able to infiuence essentially on the a-decay 
dynamics and, therefore, to change all its characteristic. From this point of view, the discussions |141 115| open 
a way for obtaining a new information about the a-decay — through the angular analysis of the bremsstrahlung 
during the a-decay. In such idea we see the first possible indication to study infiuence of nuclear deformation 
on the bremsstrahlung spectra and in this direction it can be interesting to study bremsstrahlung in a-decay of 
strongly deformed nuclei. But for such researches a model describing the bremsstrahlung in the a-decay, which 
takes into account the angle between the directions of the a-particle propagation (or tunneling) and the photon 
emission and is based on the realistic a-nucleus potential, is needed. 

In [TBI Hi] the next developments in non-stationary description of the a-decay with the accompanying 
bremsstrahlung and calculations were presented. For this time, one can refer to papers [T51 \Wi HOI [H] with 
study of dynamics of tunneling in the a-decay; an effect, opened in |22| and named as Miinchhausen effect 
which increases the barrier penetrability due to charged-particle emission during its tunneling and which can 
be interesting for further study of the photon bremsstrahlung during tunneling in the a-decay. 

In |23j we developed a multipolar quantum method with taking into account the magnetic component Ml. 
While here there was a difficulty in obtaining convergence in computer calculations, we at first time found a 
way to study the bremsstrahlung angular distribution in a-decay in fully quantum approach. It turns out that 



computer calculations of the angular bremsstrahlung spectra are essentially more complicated, if to pass from the 
a-nucleus potential used in [71 1101 123| to nuclear realistic one. In [24] we constructed a new quantum approach 
for the calculations of the angular bremsstrahlung spectra where the angle was determined on the basis of a 
simple idea (started in [25]) and where at first time the realistic nuclear component of the a-nucleus potential 
was used (in the form |26|V With a resolution of the divergence problem (which plays a key role in obtaining 
of the reliable spectra) by such approach we estimated the spectrum for ^^°Po obtaining a little difference with 
results [71 llOj , and as we find such estimation is in the best agreement (in fully quantum approach) with later 
obtained experimental data ^27| for this nucleus (see Fig. 1 in [21], Fig. 5 in [27 , Fig. 3 in _28]). 

Then in [29] a new mechanism of formation of electromagnetic radiation that accompanied a-decay and was 
associated with the emission of photons by electrons of atomic shells due to the scattering of a-particles by 
these atoms (polarization bremsstrahlung) was proposed by M. Ya. Amusia, B. A. Zon and I. Yu. Kretinin. 
It has shown that when the photon energy is no higher than the energy of K electrons of an atom, polarized 
bremsstrahlung makes a significant contribution to the bremsstrahlung in a-decay. 

With a purpose to study deeper the bremsstrahlung emitted in the a-decay of ^^"Po discussed in [HI [T5] . 
its new high-statistics measurement '27' was performed for the photon energies up to about 500 keV. Authors 
found the measured differential probability in good agreement with theoretical results obtained within the 
semiclassical approximation as well as with the exact quantum mechanical calculation. Here, it was shown that 
due to a small effective electric dipole charge of radiating system a significant interference between the electric 
dipole and quadrupole contributions occurs. With a purpose to find a unified description that incorporates both 
the radiation during the tunneling through the Coulomb wall and the finite energy E-^ of the radiated photon 
up to E^Qal ^/v (where Qa is the a-decay Q-value and u is the Sommerfeld parameter), this group proposed 
in |30j a semiclassical theory of a-decay accompanied by the bremsstrahlung with a special emphasis on study 
of ^-'^"Po. The corrections with respect to previous semiclassical investigations were found to be substantial, 
and good agreement with a full quantum mechanical treatment 7, was achieved. Here, authors found that a 
dipole-quadrupole interference significantly changes the a-7 angular correlation, and obtained good agreement 
between their theoretical predictions and experimental results. 

In [28j the bremsstrahlung emission has been measured by the a — 7 coincidence to investigate on the a- 
decay dynamics of the ^^"^Po nucleus. This experiment was performed using the ^^^Ra source and the apparatus 
with Si-detector for a-particles and NaI(Tl)-detector able to collect photons with energies up to about 1 MeV. 
The experimental data with quantum mechanical calculations are found in the good agreement between theory 
and experiment for the photon energies up to 765 keV. At present day, this result is in the best agreement 
between theory and experiment in such problem. Here, at the first time a presence of slight oscillations in the 
experimental bremsstrahlung spectrum is established. In [31' the spectrum of probability of the bremsstrahlung 
emission accompanying the a-decay of ^^^Ra by measuring the a-7 coincidences and using the model presented 
in previous study on the a— decay of ^^^Po was published. These experimental data are found to be in a 
good agreement with the quantum mechanical calculations of this group. The differences between the photon 
spectra connected with the a-decay of the ^^^Ra and ^^"^Po nuclei was explained. For two mentioned nuclei 
the bremsstrahlung emission contributions from the tunneling and external regions into the total spectrum 
are estimated, and we established the destructive interference between these contributions. A phenomenon of 
the emission of the bremsstrahlung photons during tunneling of the a-particle has been established (such a 
phenomenon has been confirmed at first time experimentally for the a-decay and theoretically with taking into 
account the realistic a- nucleus potential). 

It needs to note that arsenal of experimental data is not rich. A serious difficulty in obtaining of desirable 
accuracy in measurements lies in small values of the photon emission probability. It is not clear which nucleus 
(from available ones) should be used in order to obtain the maximal probability and performance of experiment 
will be easier. One can think that theoretical estimations of the bremsstrahlung spectra should help with 
finding of nuclei and types of decays where the emission of photons will be maximal (see [H]). In aspect of 
development of the models, a main difficulty in obtaining of reliable values of the probability of the photon 
emission is concerned with slow convergence in calculations of matrix elements, caused by slowly damping 
behavior of their integrant functions. A desire to take correctly into account nuclear component of the a- 
nucleus potential, which parameters are determined on the basis of the given nucleus, intensifies this problem. 
Overcoming of such difficulty gives enough large time of calculations, that often makes further analysis of results 
hard. This have caused a necessity to revise algorithms of calculations of wave functions and angular formalism 
of the matrix elements in the multipolar approach '23] . 

In this paper the improved multipolar model of bremsstrahlung accompanied the a-decay is presented. The 
angular formalism of calculations of the matrix elements, being enough complicated component of the model, is 
stated in details. A new definition of the angular (differential) probability of the photon emission in the a-decay 



is proposed where direction of motion of the a-particle outside (with its tunnehng inside barrier) is defined 
on the basis of angular distribution of its spacial wave function. In such approach, the model gives values of 
the angular probability of the photons emission in absolute scale, without its normalization on experimental 
data. Effectiveness of the proposed definition and accuracy of the spectra calculations of the bremsstrahlung 
spectra are analyzed in their comparison with experimental data for the ^^°Po, ^^^Po, ^^^Ra and ^^''Cm nuclei, 
and for some other nuclei predictions are performed (in absolute scale). With a purpose to find characteristics 
taking influence on the bremsstrahlung probability strongly, a dependence of the bremsstrahlung probability on 
effective charge of the decaying system is analyzed. As a natural result, as supposed type of decaying system 
the emission of proton from nucleus is studied, for which the effective charge is larger in a comparison with 
a-decay. For some proton emitters estimations of the bremsstrahlung probability are obtained (at first time, in 
fully quantum approach). Also the bremsstrahlung in fission of the ^^^Cf nucleus is analyzed (at first time, in 
fully quantum approach). 

2 Motion of electromagnetic charge inside field of nucleus 

Let's consider a particle with mass m moving (with possible tunneling) inside a field of nucleus with potential 
U(y). Hamiltonian of such a system is: 

If the particle is electrically charged then it is under action of the electromagnetic field with vector potential 
A(r,t) of such nucleus. In this case, we have (see [32], p. 187): 

iJ = -^ (p - Zeff- A(r, t)) ' - Z,s^ A* X rot A(r, t) + U{y) ^H^ + W (2) 



where 



e /_ . . _\ „ e?i . „9 e^ 



W = -Zeff -— (pA + Ap) - Zeff -^^ M X rot A + Z^g - — - A^. (3) 

2mc V / 2mc 2mc^ 

Here, Z^s is effective charge of the composite system {particle-nucleus), fi is magnetic moment (we assume that 
it is very small, ^ —> 0). Taking into account Coulomb calibration (divA = 0), and neglecting item at A^/c^, 
we obtain form of the operator W: 

W = -Z,s — Ap. (4) 

mc 

2.1 Theory of perturbations of quasistationary states 

We shall study a-decay of nucleus. The decaying nucleus we shall consider as a composite system: a-particle 
and daughter nucleus. Let ^^{t) be wave function (WF) of the quasistationary state of this system with decay 
at level Ek, Ho{t) be hamiltonian of the system. Write: 

^^'-^=Hoit)^r^it). (5) 

A solution of the non-stationary wave function (satisfying Schrodiger equation at Ek) can be written so: 

*("Ht)= J ak^^°\t)dk (6) 

where ^J^, (i) and Uk are non-stationary wave function and weight amplitude for the decay at level Ek- 

a-particle during its motion inside electro-magnetic field of the daughter nucleus emits photons. We assume 
that this process is possible also during tunneling of the a-particle. Process of the photon emission changes total 
energy of the studied system and is considered as perturbation of such system. We shall study a spontaneous 
emission of photons, i. e. when the field acts on the a-particle till the emission of the first photon. Therefore, we 
deals with the perturbation acting on the system during finite period of time. Here, time moments of beginning 
and finishing of the action of perturbation on the system are defined by such beginning and finishing of time 
interval, when we assume that the photon emission by a-particle is possible. Write the total hamiltonian with 
perturbation so: 



Hit) = Ho{t) + Wit) (7) 

where Wit) is operator of perturbation dependent on time. 

The problem lays in approximated determination of new wave functions ^(i) on the basis of wave functions 
^k (^) '-'^ quasistationary states of the unperturbed system. We shall be looking for the unknown wave functions 
solving the non-stationary Schrodinger equation with the perturbed hamiltonian 

,n^ = {Hoit) + wit))^it) (8) 



^it)= I akit)¥^\t)dk (9) 



in the following form: 

where coefficients afc(i) are already functions of time (like sec. 40 in |33j, p. 177). Substituting ^ into ([5]), we 
have: 

,(0) 



^^^^°Uk + in J akit)^dk^Ho J a,it)¥^^dk + Wit) J 0^(0*1°' 

f^j-nin kj-nin kmin k-min 



Taking into account (O, we obtain: 



,(0) cMt) ^^ _ [ TJ/^,^„,^,^,T,(o) 
dt 



^^ I ^T^^dk^ / Wit)a,it)^^^Uk. (11) 



Here, we move operator Wit) under the integral, using property of linearity (superposition) of its action on the 

set 

to: 



set of functions '^i. ' with different k. Assuming the weight amplitudes a-^^it) to be dependent on time, we come 



Multiplying both parts of this equality by ^m on the left and integrating over whole volume r, we obtain: 

, j (^j ^>t'^*^>fAr^ -^dk^ I a,it). (I v,(0).H^(i).,(0)dr) dk. (13) 



fc 



Let's assume that the decaying system is described by such quasistationary wave functions ^^ which are 
satisfied to the following condition of normalization (in continuous energy spectrum; like sec. 5 in [33j . p. BO- 
SS): 

y'vl,W,*(i)^(0)(i)dr = 5(fc-m). (14) 

We rewrite (TTSl) so: 



lh^^^= J akit)Wrnkit)dk (15) 

kmin 

where 

WMt)= J ^i°^'* Wit) ^^k^dr (16) 

are matrix elements of the non-stationary perturbation. 

As the unperturbed wave function we shall use wave function of the quasistationary j-state before the photon 
emission, which coefficients in ^ correspond to: a] — (5(fc — i) (i. e. a^ = at A: ^ z). To obtain the ffi-st 
correction, we shall be looking for a^ in the form Ok — '3;[. + a^ where we substitute ak — a^, into the right 
part of the equation (fT5|) (already having values Wmk)- This gives: 

,n^^^ = Wkdt). (17) 



To point out, which of the unperturbed functions the correction is calculated to, we add the second index to 
the coefficients 0^: 



^^{t)^ I ak^{t)^f\t)dk 



(18) 



So, we write result of integration of (|17p in the form: 

By such weight amplitude wave functions in the first correction are determined. 



(19) 



2.2 Matrix element of transition 

Let's consider wave packets (WPs) of the form: 

+ CXD 



*<,/(r,t) = J g{k-hj)^,j{k,v) e-™W* dk. 



(20) 



We shall use them as definitions for non-stationary wave functions in the initial and final states. For such 
packets we have the following properties: 



-iHot 



°*^,(r)=V^(r)e-™'^ 



(21) 



Now we shall assume that the system in the initial and final states is characterized by a number of photons 
emitted. For accurate description, we come from consideration of the system without emission, which is included 
into the hamiltonian as a perturbation by the external field, to consideration of the system with presence of 
emission where the initial and final states are defined by wave function already dependent on numbers of photons. 
As before, we define the matrix element of transition i — > / of this system as the first correction of perturbation 
(fT9| (like sec. 41-42 in [33]): 



af^{t) = -- / (fc/,7ife + l\W{Y,t')\h,nk) dt' 



(22) 



where ^i(r, t) = |fci) and ^/(r, t) = |fc/) are non-stationary wave functions in the initial i-state and final /-state 
which do not contain number of photons emitted, hw = E, Uk is number of photons of one sort with impulse 
k in the initial j-state. Operator of interaction W has the form ^, Z^s is effective charge of the system, m 
is reduced mass of the system and A is vector potential of the electro-magnetic field of the daughter nucleus. 
If rifc = than emission is named as spontaneous. If n^ > than emission is named as induced. A presence 
of photons in the initial state stimulates additional emission of other photons of the same sort (this follows 
from formula of the emission probability) . The matrix element (|22p is defined in the first approximation of the 
perturbation theory. One can develop formalism simpler in the system of units where h — 1 and c = 1, but we 
shall write constants h and c explicitly. Rewrite the matrix element a/i so: 



af.,{t) = -- /(*/(r,i'), nk + l 



W{v,t') 



*,(r,t'), nk)dt' 



+ C)0 



g{k2 - kf) V'/(fc2, r) 6-™^=^)* dfc2, Uk + 1 



to 

+ 00 +OC 



V7(r,0 



-|-oo 



.g(fci -A:,)^»(fci,r)e-™('=i)* dfci, nu) dt' - 



-n J *^ 





dki g*{k2 - kf) giki - hi) ■ I ^pf{k2,r), n^ + 1 



^iw(k2)t' ^(i.,i')e-™('=i)*' dt' 



to 



■0,(fci,r), Uk ) = 



dk2 / dki g*{k2 ~ kf) g{ki ~ h) ■ [ /, /c2, nk + I 




to 



™(fc2)t'||/(j.^^')g-»«-(fcl)t' ^^' 



i, ki, Uk 



6 



or 



where 



+OC +00 

<^fi{t) = / d.k2 / dki- g*{k2- kf)g{ki~ ki) ■ (k2,nk + l\W{r,t)\ki,nk) 




W{r,t) 



,4102* 



W^(r,i')e"™^* dt'. 



(23) 



(24) 



Now let's use the following form of the vector potential A(r,i) (see p. 22, 28 in [51]): 



A(r, i) = E ( 2k,aAk,. + c+„A*^ ) , Ak,„ = ,/M£!e(")e'(k'-™*) 



(25) 



where e^"' are unit vectors of polarization of the photon emitted, k is wave vector of the photon and w — k — k . 
Vectors e^") are perpendicular to k in Coulomb calibration. We have two independent polarizations e*^^^ and 
e*^^) for the photon with impulse k {a ~ 1,2). Taking into account ([25|) . rewrite ^ so: 



mc ^-^ \ w \ ' J 



dr 



(26) 



Substituting this expression into ((24)) . we obtain: 

t 



7, „ ^ 



The transition matrix element transforms to: 



^k.a*^ 6 



(a) i{kr~wt) , ~+ (q),* -i(kr-tot) 



c^^e' " e 



1*' ^.' . l_ 



dt 



dr 



(27) 



{k2,nk + l| M^(0 I ki,nk) = -- ( fc2,nfc + 1 



g4«'2t v{r,t')e~''"'' dt' 



ki,nk ) = 



■'eff 



2_^\l ( k2,nk + l 



Ckae^^'e*"""" / e*("'2-»i-»)t' ^^/ _^ 



, d 



dr 



Q+ g(Q),*g-Jkr / piiw2-w 



r) 

vi+w)t ../ _^ 



dr 



ki,nk 



Photons are under Bose statistics, and we have: 

(nfe + l| c^^^ I Uk) = \/nk + 1, 
(rifc + l| Cfe,Q I rife) = 0. 

According to this, wc exclude the first item from ([25)1 and obtain: 



(28) 



(29) 



(k2,nk + l\W{t)\ki,nk) ^Zeff— V 

\ I I / '' '' m ^ — ^ 



k.a 



2TTh 



■(ko 



e(a)^*g— 'ki 



d_ 
dr 



t 
fcA-V^I^rfT /"e*("'^-'"i+"')*' dt'. (30) 



We shall study the spontaneous emission, i.e. emission of one photon with impulse k, before emission of which 
there were no other photons of this sort {rik = 0). We omit summation by n^: 



{k2A\W{t)\k,,Q) = Z,fl^ J2 



a=l,2 



2TTh 



■ kn 



e(")>*e^*"" 



d_ 
dr 



t 

ki\ ■ /"e*("'^-™i+'")*' dt'. (31) 

to 



2.3 Direct indication on a presence of oscillations in matrix element in non- 
stationary approach 

We note a physically interesting result, following from the formula (|3ip . If we have the perturbation acting on 
the system (a-particle and daughter nucleus) during the finite period of time, then we have the finite upper 
limit t of time integral and can use its lower limit to = 0; ^ind obtain (in case wi — W2 ^ w)'- 



Q=l,2 



= z, 






eff 



a=l,2 



2TTh 
W 



e(a)'*g-''kr 



2TTh 

w 

d_ 
dr 



r) 
dr 



fci 



^i{w2—wx+w)t' 



i{'W2 — Wi+ w) 



t'=t 



fcl 



sin {w2 — Wi + w) t — i cos {w2 ^ wi + w)t + i 



or 



{k2,l\W{t)\ki,0) = Zeff- V \—-{k2 

a=l,2 
Zeff - E \/^ 



e(a)'*g-ikr_^ 



d_ 
dr 



W2 — Wi + W 



fcl 



m 



a=l,2 



r) 
dv 



fci 



W2 — Wl + W 

sin {w2 — wi + w) t — i cos {w2 — wi + w) t 

W2 — Wi + W 



(32) 



Here, the first item is stationary function (independent of duration t of the perturbation action) which usually 
determined completely probability of the photon emission in the a-decay in stationary approaches (which at 
such definition has monotonic behavior). But the second item is oscillating function in dependence on w and 
it includes oscillations into the total matrix element. Period of such oscillations is connected directly with 
duration t of the perturbation action. The probability defined on the basis of such total matrix element has 
both monotonic and oscillating components. In result, we theoretically have obtained (at first time) a direct 
indication on a presence of oscillations in the bremsstrahlung spectra in a-decay! In such context, period of 
these oscillations has information about duration of location of the a-particle inside the region of electromagnetic 
forces of the daughter nucleus. 



2.4 Stationary approximation 

Now we use the following limits of time integrals: 

tf) = — oo, ti = +00. 
Taking into account the property: 



+00 



e"^'dt = 2TT5{a), 



from (I3T]) we obtain: 



{k2,l\W\k„Q)=Z,ff^ E 



a=l,2 



27rh 



e(a),*g-«kr 



dr 



fci ) • 2'kS{w2 — Wi + w). 



(33) 



(34) 



(35) 



Using the following notations: 



Ffi = Zeff — \l P{h,kf), p(fc,,fc/) = ^ e(")^*p(fcj,fc/), p(fc,,fc/) = /fc2 

we rewrite expression psp so: 

(fc2,l|VF|fci,0) =F2i ■2tvS{w2-wi+w) 
and the total matrix element (I^H)) has obtained the following form: 

Ofi = / dfc2 / dfci 5*(fc2 — fc/) g{ki — fci) • F21 • 2TrS{w2 — wi +w). 



dr 



ki\ (36) 



(37) 



(38) 



For quasimonochromatic packets we have: 

af, = (Afc)2|C7|2 . Ff, ■ 2TrS{w2 - Wi + w) (39) 

where C is constant. We define it, using the foUowing norniahzation for the quasimonocroniatic packet: 

(Afc)2c2 = 1. (40) 

Then we obtain the matrix element in the following form: 

"■fi = Pfi ■ 27r 5{'Wf -Wi+w). (41) 

This expression coincides (up to factor 27r) with a general definition of the matrix element in approach of 
quantum field theory (for example, see [35] (21.2) sec. 21, p. 168-169 where function Ffi for hound diagrams is 
smooth and hasn't other (5- functions). Singular factor in (j4ip corresponds to conservation law of total energy 
of the system with emission. 

Now we shall be interesting in probability of transition defined on the basis of square of the matrix element 
Ofi. Usually, in quantum mechanics the probabilities of transitions are defined in time unit and in unit of spacial 
volume. In quantum field theory they are defined if the matrix element has 4-dimensional (5-function. In our 
case, the matrix element (|4ip has one-dimensional (5-function and direct calculation of its square does not give 
the probability in space unit. However, one can resolve this question introducing further absolute probability 
(see sec. [6]) like passing from the probability to cross-section in collisions theory. 

In calculation of square of the matrix element we deals with product of two delta- functions which are singular 
functions. Such product we calculate using approach of quantum field theory (we take it from [35] . sec. 21, 
p. 169) used for determination of the transition probability on the basis of matrix element. Taking into account 
that one-dimensional (5-function appears in result of integration by whole time interval (T —> c»), we find 
formula of power reduction of d -function: 

[6iw)f = 6{w) (5(0) ^ S{w) (27r)-i f dt ^ ^(w) {2tt)-^ T (42) 

and obtain: 

\af^\^ =2TTT\Ff,\^-Siwf-w^ + w) (43) 

that looks like (4.21) in 35j (with accuracy up to factor (27r)^, see sec. 21, p. 169). Singular factor T disappears 



after passing from ()43p to the transition probability in time unit. 

We can find square of the matrix element (|4ip by another way — in approach of pure quantum mechanics (like 
sec. 42, p. 188-189 in [33]). Let's write the matrix element before application of the stationary approximation 
so: 

T 

aj,^-'- jWfdt)dt. (44) 



We shall find Wfi{t) for quasimonochromatic packets with the spontaneous emission. From (|23p. (|3ip . (|l(7)l we 
obtain: 



m ^--' y w \ dr / 



Q = 1.2 



(45) 



where the function Ffi is defined in ([36| . Substituting W/i into (|44| . we obtain (at to — 0, ti — T): 

^i(wf-Wi+w)T _ 2 

Of, = -I Ff, (46) 

Wf ~ Wi+ w 

and for square of the matrix element we find: 

. 2 (^/ - Wi + w) T 
4 sm ^^ — 

l«/.P = \Ff^\' 7 f^^ ■ (47) 

(Wf — Wi + wy 



We use the following formula (see p. 188 in [33]): 

sin (\T „, , , , 

lim^5— = 7rJa 48 

and at large T — > +00 we obtain: 

K.I- -^T|^,.r^( -/ -'+- ). (49) 

Using formula: 

S{ax) = ^, (50) 

|a| 

we obtain: 

|a/»p = 27rT |F/,|2 • 5{wf - w^ + w). 

So, we have obtained formula (|43| which coincides with (42,5) in [33] (see p. 189) with accuracy up to factor 
Ti. However, factor fi (introduced after use of the emission operator ([26])) is included explicitly into the matrix 
element Ffi which gives coincidence between ([43]) and formalism in sec. 42 from [33; complete. 

3 Multipole method of calculation of matrix element p (/cj, kj) 

3.1 Linear and circular polarizations of the photon emitted 

Now let's consider the matrix element p {ki, kf). We have: 

d 



p{h,kf) = {k2 
and 



dr 



fci\=y'^;(r)e-'''"-|;^,(r)dr 



(51) 



(fc„ kf) = Y. e(")'*p(fc., fc/) = J2 «^^"^'* / ^/ W e"'"""!: Mr) dr. (52) 

--. 1 o ^ 1 n 'J 



P 

Q=l,2 a = l,2 



Rewrite vectors of linear polarization e^"-* through vectors of circular polarization ^^ with opposite directions 
of rotation (see [33], (2.39), p. 42): 

e-i = ^(e«-*e(2)), ^^^^_i_(e(i)+,e(2)), ^o = e(3) = 0. (53) 



For vectors ^±1 we obtain: 



1 /. . \ \ r M^. 1 



Wc find: 



e(^)^^(e..+H ; [^'^^■*-7iu;^+^- 



(54) 



E e(")- = --L(e;, -Ci) - ^(Ci +ei) = ^Ci - ^Ci = /^-iCi +/^+iai (55) 

where 

l-i , 1+i , , 

For numbers /i±i we have: 

h^i + h+i = -iV2. (57) 

Then one can rewrite p {ki, kf) so: 

p{h,kf)^ J2 h,Q frf {r)e-^'^'^^ Mr) dr. (58) 

M=-l,l 
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In further calculations of the matrix element p {ki, kf) the different expansions of function e^*'^'" connected with 
the vector potential A of the electro-magnetic field of the daughter nucleus are used. On such basis one can 
construct different approaches for calculation of the bremsstrahlung spectra. Important reason of the expansion 
of e"**^"" lays in transition from the initial integral in (j58[) . convergence of which is less possible, to series of 
higher convergent integrals. Usually, the first such integrals are used in estimation of the spectra of the photons 
emitted. 



3.2 Expansion of the vector potential A by multipoles 

We expand the vector potential A into multipoles. According to [31] (see (2.106) in p. 58), we have: 



^ue 



ikr 



/2^ ^ (2; + 1)1/2 -i jji^^^^^ Q^ 0) . r^^^^j.^ ^^) ^ ^^iAUv, e) 



(59) 



where (see (36], (2.73) in p. 49, (2.80) in p. 51) 

Aj^(r,M) = ji{kr)Tii^^{-aph) 
Ai,ir,E) = 



2^-pYj/_i(fcr) Tu-i,r.inph) - y 2rTT^'^^^'^^^ Tu+i^^{nph). 



(60) 



Here, Ai^(r,M) and Ai^(r,E) are magnetic and electric multipoles, ji{kr) are spherical Bessel functions of 
order I, T/;/_y(n) are vector spherical harmonics, 9i, 62, Os are angles defining direction of vector k relatively 
axis z in selected frame system. According to [36] (see p. 45), the functions Tw i,{ri) have the following form 
(Co = 0): 

Tj7,m(n) = ^ {l,l,j\m~ fi,fi,m)Yi^„i^^{n) ^f, (61) 

where {l,l,j\m — fi,fi, m) are Clehsh-Gordon coefficients. 

Formula (|59p is defined when the vector k is oriented arbitrary concerning the selected frame system. If to 
use such frame system where axis z is parallel to the vector k then from ([55)1 we obtain (see [3S], (2.105) p. 57): 



^^.e*"- = //y2^ ^ (2; + 1)1/2 z' rAi^(r, Af) + z/iA,^(r,S) 



(62) 



This formula is useful when it is not important to orient the frame system relatively the decaying nucleus. This 
is the case when the nucleus and decay are spherically symmetric. If the decay is asymmetric relatively axis z 
then it needs to orient the frame system relatively the decaying nucleus. Here, the vector k and axis z are not 
parallel and it needs to use ([59]) . 



3.3 a-decay in a general case 

At first, we shall consider a case when the decaying nucleus has a general form. From (|58p and (|59)) we find: 

p{h,kf) = 



J2 /iM/v^Kr)-A^^^E(-*)'^/2m-i?[.';(^,0,o); 

— 11'' / , . 



^^-1,1 



[Al{r,M)-tf,Al{r,E)] . f_^,(r)j • dr 



Here, matrix-function D\)*{lp, 9, 0) defines direction of vector k relatively axis z in the frame system for r: angles 
Lp and 9 point to direction of vector k, but not the vector r. So, we obtain: 



p(fc„fc/) = V2^^(-z)'\/2m Y. h^pi-D'f^{^,9,Q) 



^=-1,1 



Pfy -it^pt 



where 



Plv 



Ply 



Jrf{r)(-^Mr))Al{r,M)dr, 
j rf{r)(^Uv)\ Al{v,E) Av. 



(63) 



(64) 
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3.4 Approximation of the spherically symmetric a-decay 

Now we shall study a case when the decaying nucleus can be considered as spherically symmetric. In such 
spherically symmetric approximation, wave functions of the decaying system in the initial and final states are 
separated into radial and angular components: 



V'i(r) = ipiir)Yoo{nl.), 
V'/(r) = ^f{r)Yi^„,{nf). 

We rewrite integrals ([64|l so: 

+00 

drjdnr'rfir) [^^Mr^j AL(r,M), 

O 

drjdnr'rfir) (|:V'.(r)) AUr,E). 



(65) 



Pliy — 



Plu — 





+00 



(66) 



In expansion of the potential A one can use formulas ([5^ and ([SO]) . Taking into account: 

Di^{^,e,Q)^5,,, (67) 

from (l63l) we obtain: 



p{h, kf) = V2^Yl (-*)' \/2m pf - ipf (68) 

I 

where 

Using the following formula (see (2.56), p. 46 in [5S]): 



I /(.) vi„.,„,, ^ y^ji^ (I + ^/ T„.,,„.(„,-) - ;^ (I - ^/) T„«„(„.) (70) 



and taking into account (|65p . we obtain: 



^^^«--^T„,oW). (71) 



Now we shall calculate expressions ([S^ . For the magnetic component pf we obtain: 

+00 
pf = ^ /i/^M / drjdnr^rfir) {^^Mr^j A;;(r,M) = 

M^-1,1 

+ OC 

- Y. ph/j dr j dnr\){r) Y^^^M) (^-^^T^MK)) K.ikr) T,;,,^,^^(np,) 

+ OC 

= - f ^}{r)^^KAkr)r^dr ■ ^ ^^/i^ / r,;„(n,0 Toi,o(n;) T,;,^,^^,^(np„) dO. 
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For the electrical component pf we obtain: 

-foo 

+ OC 

= Y.^^^ j dr J dnr^^}ir)Y*,M){~^^Toi.o{K)y 



^^—^JipH-iikr) Tl^i^^^^Juph) - J ^^ f+i ^'p^+i^^'^) '^Ini^H+hMph) 



''ph ~r -L V P^ 

+00 



lph + ^ f :,, .d(pi (r) 



T / ^}(^)^^K.-dkr) r'dr ^ h,jY;^_^{nf)Toi,oiK)Tl^i^^_,Jn,^)dn- 



2lph , ^ . ^. ^^^^ 



^p/l 



T / 'P}i^)^^K.+iikr)r'dr ^ /i^ I y,;,„(n/)Toi,o«) Tr^^,^,^+i,^(np,) df}. 



So, the components p*^ and pf have the form: 

+CX) 



Ptl = - / '^}M^J'.'.(^'') ^''^^ E ^^h,jY*^{nf) Toi,oK) T,;,z,.,,(%.) df^, 



P; 



+OC' 



T / v}{r)'^^K,-i{kr)r^dr ^ /x^ | y,;„(n/)Toi,oK) Tr^,;^,_i,^(n„0 rff7- 



2?ph , ^ . ^. ^^^^ 



(72) 



^=±1 



Let's introduce the following symbols: 

+00 



J{lf,n) = I ^*{l,r)^^Ukr)r'dr, 



Li -^ 



(73) 



Then expressions (|72p are written so: 
PLf,^ = -lM{lf, Iph, Iph) ■ .7(1, 1), 



pf = 



-J-^f—-YdEiifJph,iph ~ 1) • JihJph - 1) + y 2; ^+ ^ iE{if,iphJph + 1) • J{ifJph + !)• 



(74) 



Using the following value of the Clebsh-Gordon coefficient (see Appendix 18.31 also ^36], Table 1 in p. 317): 

(110|1,-1,0) = (1101- 1,1,0) = yi, (75) 

from dni]) and dZID we find: 

Toi,o«) = Y. (1101 - ^^0) yi,-^(n;) Cm = y^ E ^1,-mK) C» 

/^— — 1,1 
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and from (|73p we obtain: 

/^— ±1 /l' — ±1 

If to study photon emission into all final states with different values of m and the same number Z, then as 
the wave function of the final state amidst ipf (r) in ([65|) we shall use superposition at all states with all possible 
m: 

i,f{r)=J2^fir)Yi^^inf). (78) 

m 

Let's assume that the radial component of wave function ipf{r) does not depend on m for selected If. We 
rewrite ([75]) so: 

^fir) = ^f{r) Y^YiMnt). (79) 

m 

The transition matrix element into the superposition of all states with different m at the same If is defined by 
(j74p where it needs to change expressions (j77p for the angular integrals I^^lf, Iph, n) and iMilf, Iph, n) into the 
following: 



fi=±l m •' u'=±l 



hiilfylph^n) 



lE{lf.lph,n) - J- > > /i^ >^* K) > Yi,_^,(n;)Cp'TL„ „(%/.) df] 



u=±l m -^ u'=±l 



(80) 



3.5 Vectors n^, n^, rip/, and calculations of the angular integrals 

Let's analyze a physical sense of vectors nj., n(. and np/j. At first, we consider the vector nj.. According to 
definition of tpi (r) , it determines orientation of radius- vector r from the center of frame system to point where 
this wave function describes the particle before the emission of photon. Such description of the particle has a 
probabilistic sense and is fulfilled over whole space. Here, angles 9i and ipi of nj, characterize orientation of this 
radius-vector concerning axis Oz of the frame system. 

The vector n^ determines orientation of radius- vector r from the center of the frame system to point where 
■0/(r) describes the particle after the emission of photon. Such a description of the particle is fulfilled over 
whole space and has the probabilistic sense also. Here, change of direction of motion (or tunneling) of the 
particle in result of the photon emission can be characterized by change of quantum numbers / and m in the 
angular wave function: ybo(nJ,) -^ Yiminl) (which changes the probability of appearance of this particle in 
different directions, and angular asymmetry is appeared). Angles 6f and (pf of n^ characterize orientation of 
radius-vector r relatively axis Oz of the frame system. L e. two vectors nj. and nl have similar sense and 
coincide in the angular integrals (|77p : 

nj, = n,^. (81) 

The vector iiph determines orientation of radius- vector r from the center of the frame system to point where 
wave function of photon describes its "appearance" . Using such a logic, we write: 

Uph = < = nl = Ur. (82) 

We use such frame system where axis z is parallel to vector k of the photon emission. Than dependent on r 
integrant function in the matrix element represents amplitude (its square is probability) of appearance of the 
particle at point r after emission of photon, if this photon has emitted along axis z. Then angle 6 (of vector 
iir) is the angle between direction of the particle motion (with possible tunneling) and direction of the photon 
emission. 

On the basis of ([82l) we rewrite the angular integrals (I77p so: 



(83) 



\[l E ^^M / Yi*„,{nr) E ^1,-M'K) VTr,,„,^(n.)dr!, 

^ -^ u=±l -J u'=±l 
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Writing from (|6ip vector spherical harmonic: 



Ti^;.n,M(nr) = ^ (n,l,;p;,|/x-/i",//',^)y„,^_^,,(nr)eM" (84) 



and taking into account orthogonality of vectors ^±i, we obtain 

lM{lf,lph,n) 



u— ±1 u' — ±1 



We use the following definition of the spherical functions Yim{0, f) (according to [33j, p. 119, (28,7)-(28,8); 
also Appendix 18. 2p : 



YUOM = (-l)"^^'\/^f^^^'(cos^)-e™- (86) 

where P["{cos9) is the Legendre's polynomial (see Appendix 18. ip . Rewrite the angular integral in ((85)) so: 

J l^:„(n,.)ri,-^'(n.)i;:^-^.(nr) dn = 

^ . .■^+i^-i-m'+i+m->/+u.-m'i /_-.y+„ .;+„+! / 2/ + 1 (? - |777j)[ /X / 2n + 1 (n - |m - ^1)! 

^ ' ^ ' V 47r (/ + |m|)!V8^V 47r (n+|^-^'|)! (87) 



X /" e*(-™-A''-p+M')v ^(^ . /" pl™l(cos6') Pl{cos0) P,!^"^'! (cos 61) • sinOdedip. 



The integral over yj is nonzero only in the case: 

m = —fi. 
Taking into account /i = ±1, we obtain restrictions on possible values of m and If: 

m = —ji = ±1, If > 1, 



and also 
We obtain: 



"> I/^-m'I = |m + //|. (89) 

, ^w+„_,'+i+^:i+^ ^,+„+i / 3(2^ + l)(2n+l) (1-1)1 jn - \m ^JTfjl 
-^'^^ ' V 32^ (7TI)! (n+|m + M'|)! ^ ^ 

X /" P/ (cos 61) Pi^ (cos 61) i^"+''' I (cos 6*) -sin 61 d6'(i(^. 
Let's introduce the following coefficient CT'f' : 



^m,' , w,+„+i_^'+^:i+iia . I , , . (2/; + l)(2n+l) (?/ - 1)! (n-|m + Ml)! , . 

and function/"^' (6*): 

Cf (^) = ^/,(cos0) PKcosO) Plr+^'l(cos0). (92) 
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Then we obtain the total angular integrals lui^fj^ph^ n) and /£;(//, Iph, ^), after integrating over cp: 

u'=±l n 



(93) 



(94) 



Now let's introduce the following differential expressions of these angular integrals: 

diM{if:iph,n) ^ _^/i_„^j;/+«+i y erf fr'io), 

Sin0d0 ^ Iflp^nJlfn y /' 

dlEilf;lpk,n) ^ h_^^ ^If +n+l Y erf f^f^'fQ) 

sinOde ^ hipHuJijn \ J 

ll'=±l 

and define the differential matrix elements dpf'^ and dpf dependent on the angle 6: 

_d^ ^ _ dlM{lfJpHM) . j(^ ; ) ^^,^^^^^ ^l,+lp. + l j(^i I ^^ ^ ^»V ^;V (g) 
dPl iph + ^ dlE{lf,lph,lph — i) J/, , i\ , / Iph dlE{lf,lph,lph + i) m , , i\ 



sin 0d0 y 2lph + 1 sm0d0 y j ^ f^ > ■ \i 2lph + 1 sin 61 d0 



h-rm^^^'''^^\^^^J{hJp. + l) E ^"l...+iCt+iW- 



I /i' — ±1 

(95) 
One can see that integration of such functions by angle with limits from to tt gives the total matrix elements 
pf and pf exactly. 

For transition into superposition of all possible final states with different m at the same If instead of (|95p 
we obtain: 



AiiL ,..+. 



-+ij(/,,/,.) y: ^^h^^ E cTsUjfs'>)^ 



sm0d0 

rn— ±1 fi' — ±l 



-j{ij,ip,-i) y: h-m E ^;;,c...-iCC-iW- oe) 



sin 6*^0 V 2Lh , ^ ^1,1, 

V '^ m=±l /i' = ±l 

V ^ m=±l /j'=±l 

3.6 The angular and integral matrix elements at the first values of If, Iph 

We shall find the matrix element at the first values of If and Iph- We use 

lf = l, lph = l. (97) 

From ([96]) we write: 
dpf' 



m— ±1 /i' — ±1 



dpf 



_E 



sin d0 



y|j(l,0) 5] /._„ Y. Gno'/iT'W +/^-^(l'2) E ^- E ^n2'/iT'W- 

771— ±1 /J.'— ±1 m— ±1 /J,' — ±1 



(98) 
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Calculating coefficients Cj^^ and functions /{^^ (6) (see Appendix 18. 4[ [83]) . we obtain: 
dpf' 



sinede 



= - 2. j(f , f ) • ( -h, c^,Y.fii-\e) + h^, cli-,'flf\d)j = 

3 / f 3 / 1 

= - iJ{l,l) ■ [ hi ■ - \ -— ■ sin^ 9cos9 + h-i • - \/ ;— • sin^ 9 cos ( 
8V27r 8V27r 



J{l,l)-sm^9cos9, 



dpi 



;=.E 



sin 9 d9 



= —I 



J(f , 0) I /iiCrik' fio''i&) + h-iClf„' fl^-\9) ) + 



' cos f = 



J(f , 2) I /iiCiVr' /r2''"'(^) + hiCi^^' f[,''\9) + h^iClf,' flf\9) + h^iCl^^ fl,\9) 



|.(l,0)(-.,.i/±.sin^.-._,.l/A.3in^, 
/._,.iyA.l,in^,(3eos^^-f) + /._,.l/A.3sin^ 



-vf''^'''^4v2^-("^ + '^-^'-""" 



J(f,2) 



3 

f /2 



or 



^ y^ • (/ii + ^-i) • 3 sin^ - i ^A . (/.^ + ;i_^) 1 sin2 (3 cos^ ^ - f ) 
W-- J(l,0)-sin2 6l + i- J - ■ J{1,2) ■ sin^ e (l - Ssin^ 9) 
dp{' 



sin 0d9 
dpf 



3 / f 

- \ - ■J(l,l)-sin'^9cos9, 

8 V TT 

'^ J{l,0)-sin^9 + i- J- ■J{l,2)-sin^9 (l-3sin^t 



(99) 



sinOde 8 V TT ^ ^ ' -' ' " 8 V tt 

Also we find the integral matrix elements. Integrating ([M]) by angle 9, we obtain: 

pf = --J-J{l,l)-fsin^9cosd-sin9de = --J-J{l,l)-fsin^9dsine^0, 



pf = i- J- J(1,0)- I sin^9-sind dd+ i- J- J(I,2) • /sin^6l fl - 3sin^ e*] -sin 6*^61 = 



= -i- \j- ^(1,0)- l{l~cos^9) dcos 61- i- J- J(I,2)- f (-2 + 5 cos^ 9 ~ 3 cos^ o) dcosO = 



■1 |2" r cos30 

^z-y-J(I,0). cos^-^- 



1 /l w, „x f „ . 5cos^9 3cos^9 



- .w§^(i,»)^{i4}+.iyi^(i,2)-{-2+|-|}^4yf{.(i,o)-iyi^.(i,2) 



or 



pr = 0, 



^f 



^^/!-H^'0)-ii)^-^(^'4 



(100) 
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4 Probability of photon emission with impulse k and polarization 

e(«) 

Now we define probability of photon emission when before emission (i.e. in the initial i-state) there is the given 
flux of particles with quantum numbers U — rrii = 0, and it needs to find average number of particles after 
photon emission (i.e. in the final /-state), scattering with the interesting angle relatively direction of photon 
emission and having given quantum numbers If, mj. Further, we shall consider two approaches for definition 
of such probability. 

4.1 Angular probability of photon emission in dependence on direction of impulse 
of the a-particle 

In formulation of the first definition of angular probability we assume that direction of particle motion (with 
taking into account its tunneling) before and after photon emission, direction of photon emission are defined by 
impulse pi of the particle before emission, by the impulse p/ of the particle after photon emission and by the 
impulse k of photon emitted, correspondingly. 

Probability of transition of the system in time unit from the initial i-state into final /-states with parameters, 
being in the given interval duj, with photon emission with possible impulses, in the given interval dvph, we define 
on the basis of matrix element a/i in the form (|4ip so: 

dW = -i^^ • dv, dv = dvf ■ dvph- (101) 

Taking into account Exp. (I43p for square of the matrix clement, we obtain (see [33], (42,5) sec. 42, p. 189; [51] . 
sec. 44, p. 191): 

dW ^2TT\Ff,\'^ 6{wf-Wi+w)-dv (102) 

where dv is a set of parameters, characterized photon and particle in the final /-state. In a general case, dvph 
denotes a set of parameters, described photon and having continuous region of values. If the emission of photon 
with given impulse k is considered, then: 

dk"^ w'^dwd^ph 
d-vn - (^ - (2.c)3 (1^3) 

where d^ph — d cos 6ph — sin 9ph dOph dipph, kph — w/c. Here, one assume that wave function of photon (plane 
wave) is normalized on the one photon inside unite space volume V = 1. Then dvph is number of states inside 
unit phase volume V dk^. 

Now we shall be interesting in such transition of the system, when the particle after photon emission 
propagates with impulse inside the given interval dpy. Let's define interval of states dvf, like ()103|) (also [33] 
p. 599; in contrast with [35| p. 171-175 — without devision on {2TTKf'): 



dp/ _ dp^ dpy dp^ _ p)dpf 
(27r?i)3 ~ {2T:hY ~ (^i^Kf 



d-f - TT^h^ - 'l.Zs = jh;^ ■ d^Pf (104) 



where dQ,f = dcos9f = svnO f dO f dip f , pf — |p/|. Substituting definitions for dvf and dvp^ into (jl01[) . 
integrating (|102p at dw, we obtain: 



Wf 



P/l^/^l ....n..o. .„..-.„. .„.-^'~^/ 



dW ^ ; ■: 3 dpfdVlphdVLf, Wf,=Wi-Wf = — '-. (105) 

(27r)° ft c"* n 

One can note, that in result of such one-dimensional integration by dw we obtain law of conservation of energy, 
which was not taken into account previously (like, as in quantum field theory after 4-dimcnsional integration 
one can obtain law of conservation of energy-impulse; for example, see [35] p. 171-175). Substituting Exp. ((36)) 
for Ffi into (I105p . we obtain: 



dW 



m^ (27r)4a2c3 ^^^^'^^^ '^P'"^^^'"'^^^'- i2nrn'c^ m^ 



Q = l,2 



dpf dQph dflf. 
(106) 
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In consideration of impulse interval dpf it needs to take into account, that niodul of impulse p/ in the final 
/-state is not arbitrary but is defined by law of energy conservation (|I05p . Using (as in [35], see p. 33): 






(107) 



we obtain 

Pf dpf = —hmdwph- (108) 

Taking into account this expression (its absolute values), from (|106p we define differential angular probability of 
photon emission so: 



dwph d flph d flf 

dW{ipf,0f) 
dwph d riph d cos ( 



dW{vf,0f) 



7 



dwph d ilpfi sin 9f d Of dipj 

dWiipf,0f) 
dwph d flph sin Of d Of 



(27r)4?ic3 
(27r)3?ic3 



m 



^nPf 



,("), 



Q = l,2 
Q = l,2 



p{h,kf) 



,(").* 



p{h,kf) 



(109) 



Taking into account, that the wave functions of the particle before and after photon emission have not the 
impulses Pij, but they have their wave vectors kij: 

P_U_ 

n ' 



kt.f = 



(110) 



we rewrite: 

dWiipf,Of) 


dW{(pf,0f) 


(27r)4 c3 m 

'72 2 7 

^e// e Wf^ kf 

(27rc)3 TO 


)^ e^-^^*p{k,,kf) 

Q = l,2 

> ; e(")^*p(fc„A:;) 

Q=l,2 


dwph d ilph d rtf 

dw{^f,ef) 


dwph d V,pfi sin Of d Of dipf 
dW{ipf,9f) 


dwph d riph d COS Of 


dwph d flph sin Of d Of 



(111) 



These expressions define probability of emission of photon with impulse k (with averaging by polarizations 
e^"^), where impulse of the particle after emission lays inside interval for the angle dflf or angle dOf. Note, 
that in definition of Ffi integration by whole space at r is fulfilled (i. e. averaging of wave functions of the 
particle before and after photon emission, wave function of photon is realized by whole space). 

The probability of photon emission is inversely proportional to normalized volume V, which can be used 
arbitrary. With a purpose to obtain the characteristics, which characterizes the process of emission and does not 
depend on V, it needs to divide the differential probability of emission dW on the flux j of outgoing a-particles 
in a-decay, which is inversely proportional to this volume V also. Write the differential probability so: 



dWiipf,Of) =n,v{p^)■dP, 



(112) 



where n^ is average number of particles in time unit before photon emission (in used by us normalization for 
wave function in the initial i-state we have Ui = 1), v{pi) is module of velocity of outgoing particle in the frame 
system, where colliding center is not moved (which coincides with laboratory frame, where the second particle 
is not moved). Factor P is proportional to the element of the angle of the particle after its scattering in result 
of photon emission and we shall name it as differential absolute probability (while value dW we shall name as 
the relative probability). Such a definition coincides with definition (3.4.12) for the differential cross-section in 
[39] (see sec. 3.4.4, p. 162; also sec. 4.3, p. 238-242). 

According to (9.8) from [40] (see sec. 9, p. 46), for particles with finite mass we have: 



P* = 



Ei Vi 



From here we obtain: 



E, 



Vi 



v» = 



E, 



(113) 



(114) 



that agrees with v{p) = |p|/po from ^35] (at c = 1, see sec. 21.4, p. 174). Taking into account, that wave function 
of the particle before and after emission have not impulse Pij of this particle, but have the wave vector fc^j, 
on the basis of piOp we rewrite expression for velocity so: 



E, 



(115) 
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From here we obtain equation of connection between differential relative and absolute probabilities: 



dPi^f,0,) = '^^M^dWi^^^O,) ^^ 



m v(k) 



hc^ ki 



and from here we find: 

dP{^f,6f) 



^e//^^ Wftkf 



(27r)4 c3 



m 



p{h,kf] 



dwph d^ph dflf ■ 



E,. 



hc^ ki 



Z^ff e Wph E^ kf 



or 



(27r)4;ic5 

dP{Vf,0f) 

dwph dflph dilf 

dPi^f,Of) 

dwph dflph d cos 9f 



p{ki, kf) ■ dwph dQph d^f 



Z^ff e Wph El kf 



(27r)4?ic5 



Wph E^ kf 

(27r)'^ hc^ mki 



p{h,kf) 
P{kt,kf) 



(116) 



(117) 



We also define intensity of photon emission in dependence on angles by multiplication of corresponding 
probabilities on hw (as in sec. 45 j34j ) : 



dl{(pf,t 



dwph d flph d rif 

dl{^f,0f) 
dwph d ftph d cos 6f 



Zlffe^hwl^kf 
{2ttY c^ m 

Zlff e^ ^ wl^kf 
(27rc)3 m 



1=1,2 



(118) 



We see, that so defined angular relative probabilities piip . absolute probabilities (|116p and intensities (|117p 
have real values. 



4.2 Angular probability of photon emission in dependence on direction of motion 
of the a-particle 

Now we consider another way of introduction of angular probability of bremsstrahhmg. The probability of 
transition of the system (during time unit) from the initial i-state into the final /-state, being in the given 
interval dvf, with emission of photon with possible impulses inside the given interval dvph, we define so (see 
[33], (42,5) sec. 42, p. 189; [SJ, sec. 44, p. 191): 



dW = -i^^ ■dv^2'K \Ff,\^ 5{wf - w, + w) ■ diy, 



du = dvf ■ dvph 



(119) 



where dv is values characterizing photon and particle in the final /-state. If the emission of photon wiht impulse 

k is considered then 

d^k w^dwd^ph ,.__- 

''"^^ ^J2^= (2.c)3 ' (120) 

where dflph — d cos 9ph — sinOph dOph dtpph, kph = w/c. Then integrating pi9p by dw, we obtain: 



dW^ 



w% \Ff.\ 
(27r)2 c3 



d^ph dvf, Wfi ^ Wi -Wf 



Ei-Ef 



(121) 



Now we note concerning interval dvf. In definition (|119p we use matrix element Ffi which we define as integral 
over space with possible summation by some quantum numbers of the system in the final /-state. One can 
consider such procedure as averaging by these characteristics and then Ffi does not depend on them. Therefore, 
we shall suppose that interval duf in definition (J119p takes into account only such additional characteristics 
and quantum numbers of the system in the final /-state, by which integration or summation was not fulfilled 
in definition of Ffi. 
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Substituting Exp. ([5S)) of Ffi into p2ip . we obtain: 



dVF = 



^e//^'?^ 



w/j 



27rc3 



p{h,kf) 



dflph dvf — 



27rc3 



/J-7-1Z 



a=l,2 



.("). 



P(fc*,^/) 



dQphdvf. (122) 



This expression represents probability of the photon emission with impulse k (and with averaging by polarization 
e*-")) where the integration over angles of the particle motion after the photon emission has already fulfilled. 
Such probability is averaged over all possible directions of the particle motion after emission and therefore does 
not depend on them. 

For description of the photon emission with impulse k with taking into account direction n^ of motion (or 
tunneling) of the particle after emission we introduce the probability of emission so: 

We shall define differential angular probability concerning angle 6 (and differential angular probability concerning space angle 
Vl) such a function, definite integral of which by the angle 9 with limits from to t: (definite space integral over 
angles 9 and ip) corresponds exactly to the total probability of photon emission il23\) . 



Let's consider two function: 
dW{vf,9j) 
d ^ph d rif 
dW{9f) 



dWicpf,9f) 



d Qph sin 9f d 9f dcpf 
dW{9f) 



2Tr c^ nn? 
Z^ff h e wfi 



dflf 



P{kt,kf) 



d flph d cos 9f d flph sin 9f d 9f 



2n c^ ni? sin 9f d9f 



p{h,kf) 



(123) 



Introducing differential matrix elements dependent on angle 9 and space angle SI (like definition (I96p ). we obtain: 



dw{tpf,e 



f) 



d ilph dilf 
dWiOf) 



27rc3 






p{h,kf) 



dp*{k^,kf,n 



f) 



Zljjfie^ Wf, , ., , . 
o 3 T iPiki,kf) 



dflf 
dp*{ki,kf,9f) 



h.e. 



h.c. 



(124) 



dilphd cos 9 f 271 c'^ 777.2 1^ dcos9f 

One can see that so constructed functions satisfies exactly for such definition and, therefore, they can be used 
as definitions of the angular probabilities. 

The total (integrated over angles) probability of the photon emission is: 



W 



27r C^ 777^ 



Pih,kf) 



(125) 



Intensity of emission in dependence on angles we define by multiplication of corresponding probabilities on 
fiw (like sec. 45 in [34]): 



dl{ipf,9f, Qph) 
dl{9f, riph) 



/» 



27r C^ 777^ 

v-2 t2 2 2 

27r C^ 777^ 



p{ki,kf) 
p{ki,kf) 



dp*{k,,kf,nf) 

drtf 

dp*{k,,kf,9f) 



- h.e 
h.e. 



(126) 



d cos 9f 

From (|125p and (|126p we see that so defined angular probabilities and intensities are real. 

If the probability W has dimension of mass and coincides with width F then one can define inverse value t 
to it: 

h „ „. ^^27) 



f' ^ 



W. 



According to [35j (see p. 175), in considering of transition of the system from the initial state to the final one 
the value r represents half-life of this system in the initial state, i. e. before photon emission. 

Differential absolute probability we define on the basis of the found differential relative probabilities, like this 
has done in the previous section. Using (J116p . we obtain (77^ = 1): 



dP[^f,9f) 



dW{vf,9f) _ 

77.jW(k) 



27rc3 



dW{ipf,9f) 
p{h,kf 



E,, 



72 „2 



hc^ ki 
dp*{ki,kf,Q,f] 



dfl 



f 



+ h.c. }■ dilpfi dilf ■ 



E,. 



h c? ki 



27rc5 



e' WphE^ ( dp*{ki,kf,nf) 

'p{ki,kf) 77^^^ ^+h.e. 






dQf 



dflpfi dflf 
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or 



dP{^f,0f) _ Zijfc' WphE, 



p{h,kf) 



dp*{ki,kf,Vtf) 

dnf 



+ h.e. 



dP{ipf,9f) _ ^e//e^ WphE, \^ ^ dp*{k„kf,nf) ^ ^^^ 



dflph d cos 9 



f 



2TT C^ TV? ki 



d COS 9 f 



(128) 



4.3 Multipolar approach 



Let's find the angular probability p25p in dependence on the angle 9 for the first values If — 1 and Iph — 1- 
According to (|68p . for the matrix element pi (ki, kf) we have: 



pi {h, kf) = V2^- \ {-tyV^TTi 

dpi{ki,kf) _ 



sin9d9 



2n- { {-iyV^TTT 



'pf'-^pf 


] 


= 


— V 


^ J 


1=1 


- dpf 

-sin 9 d9 


^ dpf 1 
sin 9 d9\ 


} 



= -V6^Upf'+pf), 



67r < i 



1=1 



f.E 



dpf dpi 



(129) 



sin6'd6' sin9d9 



Using the found angular electrical and magnetic components of the matrix element (j99p : 

^Cj = -\\[^-J{^A)-^^'0cos9, 
sin9d9 8 V TT 

-^^ = i\\--J{l.Q)-sin^9 + i\\--J{\,2)-sin^9(l-isinh 
sin9d9 8 V TT 8 V tt V 



and these integral components (jlOOp : 



pf 



y|.{j(l,0)-lV2.J(l,2)}, 
from (|129p we obtain: 

piih^kf) = -^yi.{J(l,0)-^^/2•J(l,2)}, 



dpi {ki,kf) . Ve 



sin 9 d9 

Now we find the relative angular probability from (I124p 
dW^^^+''\9f) 



3J(l,l)-cos6l- V2J(1,0)- J(l,2)- (l-3sin2 6') \ -sin 



(130) 



drtph dcost 



'/ 



''^ii^^^pik.,kf)'^:^'-'^Ki...^^ 



27rc3 
Z^ffhe^ 


m2 


27rc3 

.V6 r„ 


m2 



dcos9 



f 



X I 



|zy|.{j(l,0)-^V2.J(l,2)}x 
|3 J*(l,l) •cos6'-\/2 J*(1,0)- J*(l,2)- ('l-3sin2 6'') I • sin^ 6* + h.e. I 



Z^ffhe^ Wf^ ( V2 
27rc3 to2 \ 8 



J(l,0)--^/2.J(l,2) 



V2J*(1,0) + J*(l,2)- ('l-3sin2 6l) - 3 J*(l, 1) • cos0 + h.e. I 



h.e. > • sin 



Zgff iT-e w 



/« 



J(1,0)--V2.J(1,2) 



J*(1,0) + ^J*(1,2)- {l-?,sin^9\ --^ J*(l,l)-cos( 



h.e. > • sin 
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or 



d rtph d COS 9 






8 TT c'^ mri 



J(l,0)--%/2.J(l,2) 



J*(1,0) + ^J*(1,2)- fl-Ssin^ei') --^ J* (1,1) -cos ( 



(131) 



h.c. > • sin 



and the absolute angular probability from (jl28p : 



d Qph d cos Of 8 TT c^ w? ki 



J(1,0)--V2.J(1,2) 



J*(1,0) + ^J*(1,2)- (l-Ssin^ei) -^ J*(l,l)-cos( 



(132) 



h.c. > • sin 



One can see that such effect takes place: while magnetic component pf^ equals to zero, its differential part 
introduces nonzero contribution into the total probability of the photon emission. One can rewrite this expression 
so: 

dw^^+^"\ef) _ dwfHOf) .dwl^^Of) 



d Qph d cos 9f 
where the first item has a form: 



dflph dcos( 



+ A 



7 



d Qph d cos ( 



7 



(133) 



dWf\0f) 
d flph d cos Of 



'eff 



Wfz 



J(1,0)--V2.J(1,2) 



J*(1,0) + ^J*(1,2)- (l-3sin2' 

v2 



(134) 



h.e. > • sin 



and determines probability of photon emission only on the basis of electric multipole El. The second item in 
(I133P has a form: 



A 



dWf' '+''■' '(Of) 
dflph dcosO 



f 



3 Z'^fj h e^ wfi 
8 %/2 TT c3 m2 



J(l,0)--^\/2- J(l,2) • J*(l,l) + h.e. [> -sin^eicosei (135) 



and determines correction to the probabihty of the photon emission in result of inchision of the magnetic 
multipole Ml. In (jl35p one can see, that the found correction can be separated into radial and angular 
components. Therefore, in dependence on the angle it takes influence on emission equally for different energies 
of the photons emitted. From here, one can see such value of the angle, when influence of the magnetic component 
on the photon emission will be minimal or maximal: 



fiO) = sin2 6lcos6l; 
df{0) 



dO 



— sin 6* cos 6* = 2 sin cos ~ sin"^ 6* = sin 6* • ( 2 cos - sin 
dO 

sinO- (2-3sin^' 



(136) 



^2 - 3 sin^ o) 



In result, we find {a,t < < n): 



1) 01=0 

2) 6*2 = arcsin 



influence is absent, 
influence is maximal. 



(137) 



5 Nucleus— Qf-particle potential 

To describe the interaction between the a-particle and daughter nucleus we use the potential given in [25] (see 
relations (6)-(10) in the cited paper, also p8l[3T| ) in the general form 



V{r, 0, 1, Q) = vc{r, 0) + VN{r, 0, Q) + vi{r) 



(138) 
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where the Coulomb fc(r, 0), nuclear VN(r,6,Q) and centrifugal vi{r) components are 

V{A,Z,Q) 



vc{r.O) = { 2^^2 >3 ";2 3^2/ ,3. , (139) 



VN{r,6,Q) 



r-rm{e)' (140) 

1 + exp ■ 



d 

Ijl + l) 
2mr 
We define the parameters of the Coulomb and nuclear components as (see relations (14), (16)-(19) in [26j) 



"'(0 = ^S^- (141) 



V{A,Z,Q) ^ -(30.275- 0.45838 Z/A1/3 + 58.270/- 0.24244 Q), (142) 

R = i?p(l + 3.0909/i?2)+0.1243i, (143) 

Rp = 1.24^1/3 (1 + 1.646/A- 0.191/), (144) 

t = /- 0.4^/(^ + 200), (145) 

d = 0.49290, (146) 

/ = {A-2Z)/A. (147) 

According to relations (21)-(22) in [5S], we also use: 

r^{0) = 1.5268 + /?(&), R{e) = R {1 + (32Y2o{e)). (148) 

Here, A and Z are the nucleon and proton numbers of the daughter nucleus, respectively; Q is the Qa-value, 
for the a-decay, R is the radius of the daughter nucleus, V{A, Z, Q, 9) is the strength of the nuclear component; 
Tm is the effective radius of the nuclear component, d is the parameter of the diffuseness; l2o(^) is the spherical 
harmonic function of the second order, 6 is the angle between the direction of the leaving a-particle and the 
axis of the axial symmetry of the daughter nucleus; (32 is the parameter of the quadruple deformation of the 
daughter nucleus. 

5.1 Spherically symmetric a-decay 

According to [37], the deformation parameter /32 for the decaying ^^"^Po nucleus is sufficiently small that allows 
us to apply the spherically symmetric approximation for the nucleus-a-particle potential (|138p - p48p and to use 
formulas (|138p - (|14ip for the calculation of the bremsstrahlung spectrum during the a-decay of such a nucleus. 
In order to obtain the spectrum, we have to know WFs in the initial and final states. In the spherically 
symmetric approximation one can rewrite the total WFs by separating the radial and angular components: 

7 \ (149) 

VP/(r,0,0) = Rj(r)Yi^^^(6,i,) = M^ Yi^^^{0,<f>). 

We find the radial components Xij(r) numerically on the base of the given nucleus-a-particle potential. Here, 
we use the following boundary conditions: the i-state of the system before the photon emission is a pure decaying 
state, and therefore for its description we use WF for the a-decay; after the photon emission the state of the 
system is changed and it is more convenient to use WF as the scattering of the a-particle by the daughter nucleus 
for the description of the /-state. So, we impose the following boundary conditions on the radial components 

initial i-state: Xi(^ ~^ -foo) -^ G{r) -f iF{r), , , 

final /-state: Xf{r = 0) = 0, ^ ' 

where F and G are the Coulomb functions. 
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6 Calculations and analysis 

6.1 Bremsstrahlung spectra for ^^°Po, ^^^Po, ^^^Ra and ^^^Cm: comparison theory 
and experiments 

To estimate efficiency of tlie definition of the angular absolute (normalized) probability of the photon emission, 
we shall calculate the absolute spectra for the ^^"Po, ^^"^Po and ^^^Ra nuclei in such definition and the proposed 
approach where experimental data exist. 

The best result in agreement between theory and experiment we have obtained for the ^^'^Po nucleus. In 
Fig. [T] the calculated absolute probability of bremsstrahlung emission for the ^^"^Po nucleus and the newest 
experimental data in [28] for this nucleus are presented (here, there is no any normalization of the calculated 
curve relatively experimental data). In calculations, we use the second definition of the angular absolute 
probability based on direction of a-particle motion. The probability we calculate by (jl32p , using approximation 
of ? = in the calculation of the matrix element p(w,i?) (because according to our estimations the next value 
I does not give a noticeable deformation of the found bremsstrahlung spectrum). The angle 1} between the 
directions of the a-particle motion (with possible tunneling) and the photon emission is used equal to 90°. The 
nucleus-a-particle potential is defined in ()138p - (|141|) . its parameters are defined in (|142p - ()148|) . Qa-value is 
7,865 keV according to [H] (see p. 63). Radial components of wave functions of the decaying system in states 
before and after the photon emission are calculated concerning such nucleus-a-particle potential in spherically 
symmetric approximation (at /32 = 0). The boundary conditions are used in form (|150|) . For this nucleus we 
also have [41] (see p. 63): b'^^ = 100 percents (brenching ratio for a-decays that populate the given favored 
daughter state, given as a percentage of all decays), T^^^ ^ = 1,6- 10~* sec. In this figure one can see that the 
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Figure 1: The calculated absolute spectrum and experimental data [28] of the bremsstrahlung probability in a-decay of 
the ^^*Po nucleus (in calculations, we have used: Rmax ~ 2000 MeV) 

calculated spectrum for ^^"^Po by the proposed approach is in enough good agreement with the experimental 
data for this nucleus inside the region from 100 keV up to 750 keV. 

In the next Fig. [2] the calculated absolute probabilities of the bremsstrahlung in a-decay of the ^^°Po 
and ^^^Ra nuclei and experimental data in [57] and [3T] for these nuclei are presented. In calculations, the 
angular absolute probability, the nucleus-a-particle potential, its parameters, algorithms of calculations of wave 
functions and their boundary conditions are defined and calculated by the same approach as for the nucleus 
^^''Po. Here, we use: Q^-value is 5,439 keV for ^^^Po and 4,904 keV for ^^^Ra, according to [41]; the angle ■& 
between the directions of the a-particle motion and the photon emission is 90°. In figures one can see that for 
both nuclei for low energies of the photons emitted the calculated spectra are located below experimental data. 
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Figure 2: The absolute bremsstrahlung probabilities in a-decay of the spherical ^^'^Po and deformed ^^''Ra nuclei and 
experimental data in ^7 and 31 for these nuclei 

but for energies from 350 keV and higher we have obtained good agreement between theory and experiment. 
One note that for both these nuclei there is less convergence in calculations of the spectra in a comparison with 
calculation of the spectrum for ^^"^Po that can be explained by larger tunneling regions for such two nuclei. 
From here one can suppose: the tunneling region for studied nucleus is larger, the convergence in calculations 
of the bremsstrahlung spectra is obtained with larger difficulty, and study of bremsstrahlung from the tunneling 
region is more difficult. In figures one can see a little difference (tendency) between the calculated absolute 
spectrum and the relative spectrum in [241 I31j that cam be explained by different combination of integrals in 
the total matrix element in result of different expansions of wave function of the photon emitted. 

We also add the calculated absolute spectrum for ^'*'*Cm, comparing it with the high limit of errors of 
experimental data in [5l|6]. In calculations we use: Qct-value is 5,940 keV 41J and i? — 90°. From this figure we 
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Figure 3: The absolute bremsstrahlung probabilities in a-decay of the Cm nucleus and experimental data [F for this 
nucleus 

see that our calculated curve is located close to the high limit of error of experimental data. One can conclude 
that the agreement between theory and experiment is not bad and this nucleus has also sufficient experimental 
and theoretical basis for further study of bremsstrahlung processes in a-decay. This indicates to effectiveness 
of the proposed definition of the angular absolute probability of the photons emitted based on the direction 
of a-particle motion, and confirms effectiveness of the developed method of the calculations of the absolute 
spectra. 
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a-decay data 


Brcmsstrahlung probability, 1 / keV / decay 


A, 


Qa, MeV 


hT, % 


^1/2, a' sec 


100 keV 


200 keV 


300 keV 


400 keV 


500 keV 


212 


7.987 


100.0 


3.0E-2 


3.0 E-9 


8.1E-10 


2.7E-10 


9.5E-11 


3.5E-11 


218 


9.881 


100.0 


1.1 E-7 


7.5 E-9 


2.5 E-9 


1.0 E-9 


4.7E-10 


2.2E-10 


222 


8.164 


100.0 


2.8 E-3 


5.2 E-9 


1.3 E-9 


4.6E-10 


1.7E-10 


7.0E-11 


226 


6.487 


75.5 


2.5 E+3 


2.9 E-9 


5.6E-1G 


1.3E-10 


3.5E-11 


9.4E-12 


228 


5.555 


72.7 


8.3 E+7 


1.8 E-9 


2.8E-1G 


4.9E-11 


l.OE-11 


1.9E-12 



Table 1: Estimated values of the bremsstrahlung probability during a-decay of the ^^*Th nucleus and its isotopes 



6.2 Bremsstrahlung dependence on Q^ and predictions of the bremsstrahlung 
probabihty during a-decay of isotopes of Th 

In [5] it was noted about current investigations of bremsstrahlung accompanying the a-decay of the ^^^Th 
nucleus. It can be interesting on the basis of the proposed approach to estimate the absolute bremsstrahlung 
probability for this nucleus. Results of such calculations are presented in Fig. [D In calculations we use: the 
angle d between the directions of the a-particle motion (with possible tunneling) and the photon emission is 
90°, Qa-value is 5.555 keV according to [41] (see p. 63). 
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Figure 4: The predicted absolute bremsstrahlung probabilities in a-decay of the ^^^Th nucleus and its isotopes 

In [31] we explained the difference between the photon emission probabilities (both experimental and the- 
oretical results) in the a-decay of ^^^Ra and ^^^Po (at first, dependence of the bremsstrahlung probability on 
the a-particle energy was analyzed in |13p: "The difference between the two sets of data can he attributed to 
the different structure of the two nuclei, which affects the motion of the a-particle inside the barrier. The ratio 
between the two sets of data of the photon emission probability dP/dE^ is strongly characterized by the different 
a-decay energy for ^^^Po (Ea — 7.7 MeV) and ^^^Ra (Ea:—4--8 MeV) concerning the shapes of the alpha-nucleus 
harriers for these nuclei." The difference between the a-particle energies for the decaying ^^^Po and ^^^Ra 
nuclei is directly connected with different tunneling regions for these nuclei, which is directly connected with 
different contributions of the photons emission from tunneling and external regions, interference terms into the 
total spectra. And we obtained the property: The tunneling region is larger, the bremsstrahlung spectrum is 
smaller. The smaller values of the calculated total emission probability for ^^^Ra than the one for ^^"^Po can be 
explained by a consequence of the fact that outside the barrier the Coulomb field (and its derivative respect to 
r) that acts on the a-particle in the case of ^^^Ra is smaller than in the case of ^-'^"'Po because the external wide 
region results for the ^^^Po nucleus larger than for ^^^Ra and therefore the 7-emission probability for the ^'^'*Po 
nucleus is bigger. In Figl4]we have seen the demonstration of this property for isotopes of Th. In the Tabl. [T]one 
can see in detail how the bremsstrahlung probability depends on Q^-value of the nucleus for different energies 
of the photons emitted. 
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6.3 Bremsstrahlung dependence on effective charge and bremsstrahlung during 
proton emission from nucleus 

As we have seen above, Q^-value of the a-decay has strong influence on the bremsstrahlung spectra. Now we 
put a question: which are other parameters having essential influence on the bremsstrahlung spectrum? 

Let's consider Fig. 4 and Fig. 7 in [43] where the calculated 7-ray emission probabilities for the spontaneous 
fission of ^^^Ci are presented. One can see that the emission probability is changed essentially in dependence 
on mass split. From here one can suppose that the other characteristic which takes influence on the emission 
probability essentially should depend on relative numbers of mass and charge of the daughter nucleus and 
the emitting charged particle. The idea proposed in [T3] about influence of the electromagnetic charge of the 
daughter nucleus on the bremsstrahlung probability reflects this property only partially, and we see that effective 
charge of the decaying system is more directly connected with such property! Now if to consider the formula 
(|132p of the bremsstrahlung probability then one can find its direct dependence on square of the effective charge, 
i. e. we have obtained the real basis for such supposition. 

One can find that the effective charge for the a-decay is smaller (and, perhaps, essentially) in comparison 
with many other types of decays (and it is else smaller for heavier nuclei; square of effective charge equals 
to 0.18482 for a-decay of ^^*Po, and to 0.16 for a-decay of ^^°Po)! The first such example of decay, which 
has larger effective charge, we find from literature: this is proton emission from nucleus. Let's estimate the 
bremsstrahlung in this decay. To describe the interaction between the proton and daughter nucleus we use 
the proton- nucleus potential in standard optical model form given in a famous paper [32] (see formulas (5), 
(8)). Performing preliminary calculations, we neglect spin-orbit and imaginary components, add centrifugal 
component and obtain: 

V{r, I) = vc{r) + wjv(r, Rr, d^) + vi{r) (151) 

where the Coulomb vcir)^ nuclear VN{r, Rn^dfj) and centrifugal vi{r) components are 

, for r > Re, 

^c{r) = { /e2 ^ ^2 . (152) 

3 — -;::7T }■, for r < Rr, 



2Rr. 1 Rl 

VN[r,RB,,dR) ^ ^_^ , (153) 



1 + exp ■ 



dR 



Mr) = ^^. (154) 

We take radii in the form: 

RR = rRA'/^ (155) 



and define the optimum proton-nucleus standard OM parameters so: 

A^/^ + 24.01), 

(156) 



Vr{A,Z,E) = -(54.0 - 0.32 £; + 0.4 Z/^1/3^ 24.0/), 

/ = {N-Z)/A, 

TR = 1.17, 

dR = 0.75 



where A and Z are mass and proton numbers of the daughter nucleus, E is incident lab energy. In search of the 
convenient nuclei for analysis, we use Table 2 in [45] and select only 4 proton emitters which have decay from 
state 2si/2: ^^'^Ta, ^^iRg, i^Tij. and i^^Bi. 

The result of calculations of bremsstrahlung probabilities during proton decay are presented in Fig. [5l In 
calculations we use angle between photon emission and proton motion equaled to 90°. In Tabl. [2]one can see 
values of some parameters of the proton-nucleus potential for the studied nuclei. From here one can find that the 
different proton emitters have practically similar effective charges, but different essentially tunneling regions. 

We conclude: the bremsstrahlung probabilities in proton decay (from state 2si /2 j have .similar order of values 
in a comparison with the bremsstrahlung probabilities in a-decay. 
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Figure 5: The bremsstrahlung during proton decay of the ^^''Ta, ^®^Re, ^^^Ir and 
used: Rmax ~ 2000 MeV): (a) — the absolute bremsstrahlung probabilities; (b) 
studied proton emitters 
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Proton decay data 


Parameters 


Turning points 


Tunneling 
region, fm 


Eff. charge. 


Nucleus 


Qp, keV 


Orbit 




i?fl, fm 


Vr, MeV 


1-st, fm 


2-nd, fm 


}rTa83 
^flrgo 


947 
1214 
1086 
1611 


2si/2 
2si/2 
2si/2 

2si/2 


210 ms 
180 ^s 
35 ms 
3.1 ^is 


6.29 
6.3517 
6.3912 
6.6546 


-60.89 
-60.8638 
-61.0585 
-61.8599 


7.25 
7.32 
7.32 
7.56 


110.96 
88.96 
100.79 

74.28 


103.71 
81.64 
93.46 
66.72 


0.286259 
0.285328 
0.287924 
0.303988 



Table 2: Parameters of the proton decay of some proton emitters {E"^^ is square of effective charge) 

6.4 Predictions of the bremsstrahlung spectra during ternary fission 

Nuclear fission accompanied by light charged particle emission which is often called as ternary fission has been 
widely studied (see [44j and references cited therein) . Study of 7-emission during such process has been causing 
increased interest. Let's estimate the absolute probability of the photons emission during emission of a-particle 
from the ^^^Cf nucleus (which is the popular nucleus used in study of such problem). In calculations we use: 
parameters of the nucleus-a-particle potential are used according the proposed approach above, Qa-value is 
6.257 keV, R = 7.684 fm; Jf^^ = 81.6 percents (brenching ratio for a-decays that populate the given favored 
daughter state, given as a percentage of all decays), T'^?^ ^ = 1.0 • 10+^ sec. In the next Fig. [6] the absolute 
bremsstrahlung probability during emission of a-decay from ^^^Cf is presented. 

7 Conclusion 

In this paper the improved multipolar model of bremsstrahlung accompanied the a-decay is presented. The 
angular formalism of calculations of the matrix elements is stated in details. A new definition of the angular 
(differential) absolute probability of the photon emission (i. e. without normalization on experimental data) in 
the a-decay is proposed where direction of motion of the a-particle outside (with its tunneling inside barrier) is 
defined on the basis of angular distribution of its spacial wave function. Effectiveness of the proposed definition, 
the developed formalism of the model and accuracy of the calculations of the bremsstrahlung spectra are 
analyzed in their comparison with experimental data for the ^^"Po, ^^^Po, ^^®Ra and ^'^''Cm nuclei. Here, note 
the following. 

• The best result have been obtained in agreement between the calculated absolute probability of the 
bremsstrahlung emission for the ^^''Po nucleus and the newest experimental data in [28] for this nucleus 
inside the region of photons energies from 100 keV up to 750 keV (see Fig. [H Qq = 7.865 keV, the angle 
•d between the directions of the a-particlc motion and the photon emission is used 90°). 
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Figure 6: The absolute bremsstrahlung probabilities in a-decay of the ^^^Ci nucleus (in calculations, we have used: 
i?max = 5000 MeV) 

• The calculated absolute probabilities of the bremsstrahlung emission in a-decay of the ■^^"Po and ^^^Ra 
nuclei for low energies of the photons emitted are located below experimental data ^7\ and [3T] , but for 
energies from 350 keV and higher we have obtained good agreement between our model and experiment 
(see Fig.m Qc, = 5.439 keV for ^lOpo and Qa = 4.904 keV for 226Ra, -& = 90°). 

• The calculated absolute probability of the bremsstrahlung emission for the ^'^''Cm nucleus is located close 
to the high limit of error of experimental data [SJ [B] (see Fig. [31 Qa — 5.940 keV and ^ = 90°). 

Analyzing the formalism of the model, we establish: 

• A presence of oscillations of the bremsstrahlung probability in a-decay defined on the basis of non- 
stationary quantum approach (theoretically, at first time in fully quantum approach). In such context, 
period of these oscillations has information about duration of presence of the a-particle inside the region 
of electromagnetic forces of the daughter nucleus (from the a-particle formation inside the nucleus up to 
its passing through electrons shells). 

• A property: for any studied nucleus the tunneling region is larger, the convergence in calculations of 
the bremsstrahlung spectrum is weaker, and study of bremsstrahlung from the tunneling region is more 
difficult. 

• A principal difference between two experiments [F; "6] and [3] (opened at first time): in experimental 
data [5] [6] for the ^^''Po nucleus the contribution of the magnetic component Ml into the total spectrum 
is close to maximal while in the experimental data [3] for the ^^^Po and ^^^Ra nuclei we have obtained 
zero contribution of this magnetic component. 

• The non small dependence of the bremsstrahlung probability on the effective charge of the decaying 
system. 

On the basis of such model the bremsstrahlung probabilities in the a-decay of ^^^Th and some its iso- 
topes, the bremsstrahlung probabilities during proton emission from the '^^'^Ta, ^^^Re, -'^^^Ir and ^^^Bi nuclei 
are predicted (in absolute scale, at first time in the fully quantum approach). According to analysis, the 
bremsstrahlung probabilities in the proton decay (from Si/2 state) have similar order of values in comparison 
with the bremsstrahlung probabilities in the a-decay. On such a basis one can hope that experimental study of 
bremsstrahlung in the proton decay can be interesting. 
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8 Appendix 

8.1 Legandre's polynomials 

Legandre's polynomials Pi {cos 9) and associated Legandre's polynomials P™(cos0) are defined as in [33j (see 
p. 752-754, (c,l)-(c,4); [SH] (2.6), p. 34): 

Pr (cosO) = sin"' g ,/'^ ° = T^ sin" 9 ,^ ,,,^ (cos^ ^ - 1)' = (157) 

' ^ ^ (dcos6')™ 2'Z! (dcos6')'+™^ ' ^ ' 

^ ' {l-m)\2n\ (dcosfl)'-™^ ' 

where m = Q . . .1. For these polynomials the normalization condition is fulfilled (/^ = cos 6*; see [S^, (c,6)-(c,9) 
p. 753-754): 



j[Pi{l^)] 



1 

'"^^"2^' J Piitx)Pi,ifi)dfi = 0, 

-1 -1 

1 1 

|[pr(M)]'dM-^|^, Jpr{f^)pr{f^)dt. = o. 
-1 -1 

At the first coefficients I and m polynomials are: 

Fo"(cos6') == 1, 

P^{cos9) = 003 6*, 

P^{cos9) = sin 9, 



(158) 



F2°(cos^) = i(3cos2 0-l), ^^^^^ 

P2^(cos6') = 3 sin 6' cos 6*, 
Pi{cos9) = 5 81-0^9. 



8.2 Spherical functions Yi 



Im 



On the score of different definitions of the spherical functions in literature, we present the definition of them 
used in this paper. We define the spherical functions Yim{9, ip) according to [33] (see p. 119, (28,7)-(28,8)): 



Yi^ie,^) = (-l)"^^'^^f^^^'(cos^)-e™^ (160) 

where P™(cos0) are associated Legandre's polynomial (see Appendix 18. ip . For the functions Yim{9,(p) the 
following condition of normalization is fulfilled (see [33], (28,3), p. 118): 

l^r,„,(0, if) Yi^mie, ^) sin^ d9 dip = 6iv6„un'- (161) 
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ija^jlmanibm) 



TO(, = 1 



-T72- 



nib = -1 



-TJT 



1/2 



1/2 



i = ja + 1 



J =Ja 



j = ja-l 



ija + ^J^) (jg + ?•» + 1) 

(2ja + 2) (2ja + 2) y 

/ {ja + to) (jg - m + 1) 
i^ 2j, Oa + 1) 

iJa -m) (ja-TO+l) \ 
2ja {2ja + 1) ] 



ija-m){ja -171+ I) 



1/2 



1/2 



(2j 


, + 1) (2j„ + 2) 


Oa- 


m){ja +TO+ 1) 




2ja (.?a + 1) 


Oa + 


"^ + 1) (ia + m) 



2ja {2ja + 1) 



Table 3: Clebsch-Gordan coefficients. 



The functions YJmC^, <p), differed by sign m, are connected by relation (see (28,9), p. 119 in [33]): 
Write expressions for some first normalized spherical functions Yim{S,(p) (see [33], p. 754-755): 



(162) 



^00 

Yi,±i 

^2, ±1 



1 



47r 



± 



' — sin6'-e=^''^ 
_87r 

'f5 



Yw 
Y20 



I \l — cosf, 

47r 



-— (f-3cos2( 

fOTT 



COS f sm a ■ e 



■e^^f. Y, 



2. ±2 



15 
32^ 



sin2 6'-e±2^¥'. 



(163) 



^30 
i"3,±l 
i"3,±2 
y3,±3 



■i ^| cos6'(5cos2 6i-3), 

±i ^|-^ sin6'(5cos2 6l-l)-e±*'^, 

105 
i\l—-- cos 6 SIX? 6 -e^"^'^, 
327r 

±j jiL sin3 6l-e±3'¥'. 
647r 



(164) 



8.3 Clebsch-Gordan coefficients 

We define Clebsch-Gordan coefficients, according to Table .1 in [36] (see p. 317), rewriting it as Table 1. Using 
the table 121 we find: 



(011|2,-1, 1) 
(011| 0,1,1) 
(011|Q,-1,-1) = 
(0111-2,1,-1) = 



' {ia - m) {ja - TO + 1) 
(2ja + 1) (2ja + 2) 

^ {ja+m){ja +TO+ 1) 
{2ja+2){2ja + 2) 

' {ja -m){ja -TO+ 1) 

(2jg + 1) (2jg + 2) 

' {ja + to) {ja + TO + 1) 
{2ja + 2) {2ja + 2) 



' (0-1) (0-1 + 1) 
(2 -0 + 1) (2 -0 + 2) 

' (0 + 1) (0+1 + 1) 
(2 -0 + 2) (2 -0 + 2) 

' (0 + 1) (0+1 + 1) 
(2 -0 + 1) (2 -0 + 2) 

' (0-1) (0-1 + 1) 
(2 -0 + 2) (2 -0 + 2) 



= 0, 



2' 
2' 



(165) 
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(111|2,-1, 1) 
(111| 0,1,1) 
(111|0,-1,-1) = 
(lll|-2, 1,-1) = -. 

(21l|2,-l, 1) 
(211| 0,1,1) 
(211| 0,-1,-1) = 
(21l|-2, 1,-1) = 

8.4 Coefficients C,T' „ 



We define coefficients C,™f so: 



Kia-m){ja+m+l) ^ / (I -1) (1 + 1 + 1) 

2ia (ia + 1) V 2 • 1 • (1 + 1) 



'(ja + m)(ja-m+l) (1 + 1) (1-1 + 1) /T 



2ja(ja + l) V 2-1 .(1 + 1) V2' 



'(.?a-m)(j„ + m+l) (1 + 1)(1-1 + 1) /T 



2ja(ja + l) V 2-1 -(1 + 1) V2' 



(166) 



' (ja+H(ia-"^+l) ^ _ / (I- 1) (1 + 1 + 1) 
2.?a (.?a + 1) V 2 • 1 • (1 + 1) 



' (ja+m + l)(ja + m) ^ / (2 + 1 + 1) (2 + 1) ^ /3 
2ja(2ia + l) V 2 -2 -(2 -2 + 1) Vs' 



'(ja-m)(ja-m + l) _ (2-1) (2- 1 + 1) _ [l 



2ja{2ja + l) V 2 -2 • (2-2 + 1) V 10' 



'(ja+m + l)(ja + m) _ (2- 1 + 1) (2-1) _ /T 



2j,(2j, + l) V 2 -2 •(2-2 + 1) V 10' 



(167) 



'{]a-m){ja-m + l) (2 + 1) (2 + 1 + 1) /3 



2ja {2ja + 1) V 2 - 2 - (2 - 2 + 1) V 5 



^'/'--^""^ ^^ (",i,'pM "^ M,/^, mjy ^2^ (Z^ + l)! (n+|m + /i'|)! 

(168) 
Let's consider case of If = 1, ^p/i = 1 and n = Q. From ([55] ) we obtain: 

m=-/i' = ±l. (169) 



cz - (-i)^-+--'-(oii|o,;/,-™)J(^:l±lli^:il±ll(i^(^ 

"° ^ '^ V I ,/ , 7 y 32^ (1 + 1)! (0 + 0)! 



^-(01l|0,M',-m). 
Taking into account values (|165p for the following Clebsch-Gordan coefficients: 

T 



(170) 



(011| 0,1,1) = (011| 0,-1,-1) - ^1 ^, 



we obtain: 


^110 


= 0, 


1 — T^ — 






f-1-ii 

^110 


= -\ 


/ 
/ 647r 


■ (Oil 0,1,1 




"-"no 
r<ii 


= 0. 


/ 647r 


■(Oil 0,-1, 



1 / 3 



128 7r 8 V 27r 

3 1 rr 



128 TT 8 V 27r 



(171) 
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At // = 1, Iph = 1 and 71 = 1 property (|169p is fulfilled also. The coefficients C™!* ^ obtain the form: 

Taking into account values in (|f66p for the following Clebsch-Gordan coefficients: 

(111|0,1,1) = -J^, (111| 0,-1,-1) -y'i, 



we ffiid: 



C^A-' - 0, 



/ 9 ,...,„ _, 3/1 



Ci\\ = 0. 
At If — 1, Ipfi = I and n — 2 the property (|169p is not fulfilled. We have: 



(173) 



C^,' , ^.i+2+i-,'+^4^ .2111 m a' a' m) / (^ ' 1 + 1) (2 ■ 2 + 1) (1-1)! (2-|m + /x-|)! 

^"2 " ^ ^^ (211|-m-^,^,-m)W 32^ (1 + 1)! (2 + |m + ^'1)! 



y 64 TT (2 + \m + fj,'\j\ ' 

Rewrite at different m — ±1 and /^' = ±1: 



(174) 



C-i'i - f 1)1+^^ /15 (2-1-1-11)! 

^112 - (-1) y 64^ (2+1-1-11)! ^^^^1^ + ^'"^'^^ 

^^ i . (2111 2,-1,1) = 4 \/J- (211| 2,-1,1), 



64 7r 4! ' I ' ' ' 16 V 2 7r 



c^.V - (-D-^+^Ji^S^-'-^ + ^'ii .(2111 1-1, 1,1) 

"2 '^ ^ y 64 7r (2 + |-l + l|)! ^1 ' ' ^ 

^'^ (21l| 011) = -ij— -(2111011), 



64 TT ' I ' 8 V TT 



C'-' ( 111+^ / 15 (2-11-11)! I 

^"2 " ^"^^ V 64 TT (2+11-11)! (2111-1 + 1,-1,-1) 



(175) 



^•(211| 0,-1,-1) = -1^1^. (211| 0,-1,-1), 
on. - (-l)-^./^|-^-(21l|-l-l,l,-l)^ 



64 TT (2+|l + l 



^^ ^ (2111 -2,1,-1) = 1JA.(211| -2,1,-1). 



64 7r 4! ' ' ' 16 V 2 7r 

Using the found values in (|167p for Clebsch-Gordan coefficients: 

(211|2,-1, l) = y|, (211|0, 1, l) = yi, 

(211| 0, -1, -1) = W^, (211| - 2, 1, -1) = j|. 
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we obtain: 



^-11 

'-'112 



^1-1 
'-'112 



1 



16 V 27r 

1 As 



(21l| Oil) 



./_. (211 0,-1,-1) 

8 V 71^ 



8 V TT 

1 /15 



1 



1 



Ciir^ = tTfVtTT- (211 2,-1,1) = -rJ:^ 



16 V 2 7r 

1 /l5 

sV V 

1 /15 

sV V 
1 Hb 



C}h = —\ (211-2,1,-1) = —\ — 

"2 16 V 27r ^ I ' ' ' 16 V 27r 

Write the values for the calculated coefficients: 



3 

r 

io 

1 

10 
3 



1 /A 

16 V 27r' 

1 [j; 

Sy 2tt' 

8V27r' 

1/A 
16V 2tt' 



(176) 



^-1-1 

'-'110 



0, 



c^A-' = 0, 



^-11 

'-'110 



^1-1 

I' '-'110 



^-11 _ _'J / -"^ ^1-1 

^1" ^ 8'V27r' ^"1 



'-'110 — "' 



^JA C" -0- 

8 V27r' ^111-"' 



^-1-1 _ 1 / 3 ^-11 _ / "^ ^1-1 _ / "^ r*!! — / 



(177) 



8.5 Functions /™^'(6') 
We define function /™^ (0) so: 



C^'(^) = PiA^o^^) A'(cos0) /^'"+^'l(cos0). 



At Zf = 1 and n = 0, 1, 2 we obtain: 



(178) 



/i7' (6) = Pi (cos 0) Pi (cos 0) ^1"+^' I (cos 9) , 
/{7' (0) = Pi (cos 0) P^i (cos 9) Pj"+^' I (cos 9) , 



/ir (^) 



Pi {cos 9) Pi {cos 9) P; 



|m+/^'| 



(cos 6*). 



(179) 



We rewrite at different values of ?Ti = ±1 and /i' = ±1: 



fi,'^-\0) 



fw'\0) 



(9) ^ 



.1-1 

J 10 

fn\d) = 

fim 

fl2\0) = 



Pl {cos I 
Pi {cost 
Pl{cost 
Pllcost 

Pl {cos I 
Pllcost 
Pllcost 
Pllcost 

Pl{cost 
Pl{cost 
Pl (cos ( 
Pllcost 



Pl{cost 
Pllcost 
Pl{cost 
Pllcost 

Pl{cost 
Pllcost 
Pllcost 
Pllcost 

Pl{cost 
Pl (cos t 
Pl{cost 
Pllcost 



Pg{cost 
pllcost 
[cost 
[cost 



P§{c 



PoHc 

Pl{cost 
pf{cost 
pf{cost 
pf{cost 

P^{cost 
pf{cost 
pf{cost 
pf{cost 



sinf? 


sin0 


1 = sin^ 9, 


sin0 


sin0 


1 = sin 9, 


0; 






0, 






sin0 


sin0 


cos 9 = sin^ 9 cos 9, 


sin0 


sin0 


cos 9 = sin^ 9 cos 9, 


0; 







sm9-sm9-'Ss[B^9 = 'Ss[n^9, 
sin0-sin0- i(3cos2 6'- 1) = i sin^ (3 cos^ i 
sin6i-sin0- i(3cos2 6'- 1) = i sin^ (3 cos^ i 
sin 9 ■ sin 9 ■ 3 sin^ 9 = 3 sin"' 9. 



-1) 



(180) 
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